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We investigate the influence of spatially inhomogeneous chiral symmetry-breaking condensates in 
a magnetic field background on the equation of state for compact stellar objects. After building 
a hybrid star composed of nuclear and quark matter using the Maxwell construction, we find, by 
solving the Tolman-Oppenheimer-Volkoff equations for stellar equilibrium, that our equation of state 
supports stars with masses around 2 Mg for values of the magnetic field that are in accordance 
with those inferred from magnetar data. The inclusion of a weak vector interaction term in the 
quark part allows one to reach 2 solar masses for relatively small central magnetic fields, making 
this composition a viable possibility for describing the internal degrees of freedom of this class of 
astrophysical objects. 
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I. INTRODUCTION 


The study of the properties of hadronic matter under extreme conditions is perhaps the most challenging task in con¬ 
temporary nuclear physics. Due to its non-perturbative behavior, quantum chromodynamics (QCD), the fundamental 
theory of strong interactions, cannot be solved using conventional field-theoretical methods, making any prediction of 
its properties an extremely challenging task. While ab-initio lattice calculations and current heavy-ion experiments 
at the RHIC and LHC shed some light on the properties of strongly interacting matter at high temperatures, the 
opposite region of the QCD phase diagram, associated to high density and low temperature conditions, is still largely 
unknown. On the theoretical side, lattice calculations at nonzero chemical potentials are hindered by the sign prob¬ 
lem and most of the current predictions rely on phenomenological models, while experimentally none of the current 
heavy-ion colliders can reach sufficient densities to probe this region. 

In spite of the current uncertainties, the phase structure of QCD at finite densities is nevertheless expected to be 
very rich (see e.g. T]). In particular, a growing consensus has been recently building around the idea that crystalline 
phases might appear in the intermediate density region (up to a few times nuclear matter density), before the onset 
of color superconductivity. The formation of such phases could in principle delay the restoration of chiral symmetry, 
dramatically altering the properties of cold quark matter (for a recent review, see [2]). In particular, it has been 
suggested that the presence of strong background magnetic fields, a natural element in astrophysical scenarios, might 
significantly enhance the window for inhomogeneous phases Hi, possibly leading to significant effects in the equation 
of state (EoS) of dense quark matter, as we discuss below. 

While waiting for the next generation of heavy-ion colliders such as the FAIR experiment in Darmstadt and NICA 
in Dubna, which promise to access experimentally this window of the QCD phase diagram, the best laboratory for 
investigating properties of dense matter is given by compact stellar objects, which provide the only known realization 
of ultradense systems in nature. Of particular interest are measurements of masses and radii, which indirectly provide 
numerous hints on the possible EoS for QCD at high densities. In particular, the recent discoveries of stars with 
masses close to 2M© (M© being the solar mass) [5} '6] impose rather strong limitations on the possible EoS. It is 
worth recalling that a lot of modeling of the microscopic physics inside these compact stellar objects is involved when 
making this kind of prediction, introducing a large number of uncertainties. Nevertheless, while the role of hyperons 
as a softening ingredient of the nuclear EoS is still under debate IZHS1, as well as the presence of quark matter maun, 
numerous studies seem to indicate that most of the ordinary phases of (confined or deconfined) matter might not be 
able to support the large stellar masses observed. All these considerations suggest that some fundamental aspect in 
the physics of dense matter might still be missing from these calculations. Due to the high densities reached in the 
core of compact stellar objects, it might be reasonable to expect a transition to more exotic phases, whose EoS could 
be stiff enough to sustain these massive stars. Of course, the issue of the maximum mass is subject to other effects as 
well, such as high rotation rates (see e.g. [T3] and references therein) or the existence of strong magnetic fields [T51 - [TB] 
that affect the EoS and may allow those objects to support higher masses than a static, nonmagnetized star would. 

The main purpose of this work is to investigate the effects of the formation of inhomogeneous chiral-symmetry¬ 
breaking condensates in a magnetic field background on the EoS of cold and dense matter, and whether they can lead 
to predictions for compact stellar objects which are compatible with current experimental observations. In particular, 
we aim at building hybrid stars with a crystalline quark matter core, using the resulting EoS as input for the Tolman- 
Oppenheimer-Volkoff (TOV) equations. In order to build a realistic description of matter for astrophysical scenarios, 
the models under consideration will include the effects of strong magnetic fields, which are naturally expected to be 
present in the compact stellar medium. 

This work is structured as follows: in Secs. II-III, we introduce the phenomenological models employed to describe 
quark and nuclear matter. In Sec. IV, we consider a medium-dependent magnetic field and give its density profile 
for various parametrizations. In Sec. V, we build the EoS for a hybrid star with a crystalline quark matter core in 
the presence of such magnetic field, and obtain the corresponding mass-radius (M-R) plots. Finally in Sec. VI we 
summarize the results and give our concluding remarks. 


II. MODELS OF NEUTRAL MAGNETIZED QUARK MATTER WITH INHOMOGENEOUS 

CONDENSATES 


In this section we introduce the models we are going to use to investigate quark matter with spatially inhomogeneous 
chiral condensates in a magnetic field background. Since our ultimate goal is to investigate the structure of compact 
stars, we shall impose the physical conditions of electrical neutrality and /3—equilibrium in all our calculations. Vector 
interactions will also be included. For our calculations, we will consider both two- and three-flavor models, which will 
be described in the following. 
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A. Two-flavor model 


To study two-flavor quark matter in a magnetic field background we consider the following Nambu-Jona-Lasinio 
(NJL)-type Lagrangian density 

£( 2 /) _ ^ + m7° - m q ) i’ + 'ipe (i^D^ - m e ) ip e + C int , (1) 

containing a doublet of quark fields i /’ T = (ipui^d) in flavor space, with current mass matrix mn q = diag (m u ,md) 
and an electron field ip e of mass m e . A nonzero baryon density has been introduced via the quark baryon chemical 
potential fx. The covariant derivative describing the coupling of matter with a static external magnetic field along 
the ^-direction is D^ = d^ + iQA ext , with electric charge matrix in flavor space Q = diag(e„,ed) = diag(|e, — |e), e 
being the unit electric charge and A ext being the external electromagnetic four-potential taken in the Landau gauge 
A ext = (0,0,Hx,0). 

The quark interaction Lagrangian Ci n t is given by 

C-int = C\ + C 2 + Cy > (2) 

with 

Ci = Gi [{'ipip) 2 + (V5ry 5 V0 2 + W’t'V) 2 + (V’ry 5 T a V’) 2 ] , (3) 


c 2 = g 2 [WO 2 - (#7 5 V0 2 - ii’T a ip) 2 + (iprfT a ii >) 2 ], 

and vector channel 

C v = -G v [WW 2 + Wf’TW 2 + (<i’7/“>/’) 2 + W 5 7 m t£ V) 2 ] ■ 


(4) 

(5) 


Here, the matrices r a are the generators of the SU(2) flavor group and a sum in the color index is assumed in all the 
quark terms. 

For applications to compact stellar objects, one needs to consider electrically neutral and /3-equilibrated matter. 
To incorporate neutrality, we insert an electric charge term —/x e Q e in Q, with charge operator 


Qe = g VV7(W - ^dlO^d - ^elo^e- 


( 6 ) 


The electric chemical potential fx e is not an independent parameter; it has to be determined self-consistently from the 
condition of electrical neutrality. A nonzero fx e gives rise to an isospin asymmetry between the quarks, which now 
have chemical potentials 


Mm — 1 ^ ^Me > Md — M T g Me • (7) 

At this point we perform the standard mean-field approximation and introduce an inhomogeneous ansatz for 
the expectation values (ipfipf) and (ipfi"/ 5 ipf), f = u,d (in the following we neglect flavor off-diagonal mean fields 
corresponding to charged pion condensation). We note that the interaction term which corresponds to the 
instanton contribution G 2 [det V>(1+ 75 ) 1/7 + det 1/7(1 — 75 )^] jTTjUEj and is added to the Lagrangian to include the U(1)a 
anomaly of QCD, contains flavor mixing terms like (ip u 4>u' l l , d'4>d) that would significantly complicate our calculation 
when dealing with asymmetric inhomogeneous matter. As a first step, we then choose to avoid dealing with mixing 
terms by considering a simpler version of our model in which different quark flavors can be completely decoupled by 
neglecting the instanton term in the interaction Lagrangian (i.e. taking G 2 = 0). We expect the main features of the 
quark EoS to be qualitatively unaffected by this simplification, and note that at any rate the quark condensates will 
still influence each other through the neutrality condition. From now on, we will consider G 1 = Gs and work in the 
chiral limit m u = rrid = 0. 

Let us consider the following plane-wave ansatz 


-AGsiiputpu) = cos (q u z) , -4Gs(4> u i"/5ipu) = A„ sin (q u z ), 

-4 G s (ipdi>d} = A d cos (q d z) , -^Gs^dilsipd) = A d sin W). 


( 8 ) 


In the isospin-symmetric case, this ansatz reduces to the so-called “chiral density wave” (CDW) [TS], characterized 
by plane-wave condensates with magnitude A = A u = A d for each flavor and equal and opposite wave vectors, 
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q = Qu = ~qd- However, in isospin-asymmetric matter, there is no reason why different flavors should have the same 
condensate, and one should allow for two separate amplitudes A f and modulations qj, f = u,d. 

In an external magnetic field, the only rotational symmetries that survive are the SO(2) group of rotations about the 
field direction, making that direction special. That is why we have chosen the modulation of the condensates along 
the magnetic field direction. We note that a modulation along the field direction is known to be energetically favored 
in the isospin-symmetric case [19j . 

Since we are working at nonzero baryon density in a theory with vector interactions, we also introduce expectation 
values of the individual quark number densities for each flavor, 

(V’mToV’u) = Pu, {4’d'Yoi’d) = Pd , (9) 


connected to the baryon density through 3 ps = ( p u + Pd,)- While for an arbitrary spatial dependence of the chi¬ 
ral condensate a proper self-consistent inclusion of the expectation values of the quark number densities could be 
challenging, as they could be themselves inhomogeneous, in the case of a CDW modulation the procedure is actually 
straightforward, thanks to the fact that, for this particular ansatz, the quark number density is spatially constant |20j . 
The net effect of including vector interactions then amounts, just like for homogeneous matter, to the introduction of 
a shifted chemical potential for each flavor, given by mm 


Pu pu 4 Gy p u , pd Pd AGy Pd • 


( 10 ) 


Expanding around the expectation values introduced in Eqs. (|8|) and (10), we obtain the mean-field Lagrangian 


£ { mf = i’u (i'fD'p + pul 0 - A u e^ z 


V’e 


+ V’e ( ii^Dft + p^ - m, 


;) V’e 


) ipu + ipd (ii^Dj?) + Pdl° ~ A d e ll5qdZ ^j i/j d 
(Pu - Pu) 2 , (Pd - Pd) 2 


A 2 A 2 

u _ d 

8 Gs 8 Gs 


+ 


8 G v 


8 G v 


( 11 ) 


Since this Lagrangian is now bilinear in the matter fields, the corresponding thermodynamic potential can be readily 
obtained. In the following, we neglect thermal effects and work at zero temperature, a reasonable approximation 
when describing cold and dense stellar matter. The zero-temperature, mean-field thermodynamic potential of the 
two-flavor theory (11) is thus given by 


n (2/) = n e + N c ^2 n f+ 

f=u,d f=u,d _ 


8 G s 


(Pf ~ Pf) 2 

8 G v 


( 12 ) 


where N c = 3 is the number of colors and the quark contributions for each flavor are 


n f = n v f ac + ny ed , 


n v f ac = 


1 \e f H\ 


dp3 


h/A 


ds 

2 S37 2 


e~ aE f’0 + ^ e~ sE A q 

e n>0,C,e 


Ay/n (27r) 2 

V? ed = J ^Pfbf - ^ J°° dp 3 J2(\Ef,o - M/I - \E f , 0 


\reg 


|e/H| 

47T 2 


/ OO 

dp 3 'y ] ( Pf ~ Ef >n )@(pf — Ef >n ) |e=l, 

-OO __ /■ 


n> OX 


with quark energies given by 

E f, o = eyjA 2 f +p 2 3 + b f , 

Ef, n = e \j ((\f^}+P 3 + b f) +2\e f H\n, e = ±,C = ±,^>0, 


e = ± , n = 0, 


(13) 

(14) 


(15) 


(16) 

(17) 


where bf = and n = 0,1, 2,.. denotes the Landau levels. Notice that this spectrum exhibits a drastic distinction 
between the modes of the lowest Landau level (LLL), n = 0, and the rest, n > 0. The spectrum is asymmetric 
about zero for the LLL, while it is symmetric for any n > 0. The index £ is connected to the spin projection, while 
e labels particle/antiparticle energies for all n > 0. This last interpretation is not valid however for the LLL, due to 
the spectral asymmetry. Note that for n > 0, the presence of the modulation, bf 0, breaks the spin degeneracy, 
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FIG. 1. Amplitudes A/ (/ £ ( u , d )) and wave numbers bf for the CDW modulation (Eq. [8| as a function of chemical potential 
in the presence of a constant external magnetic field. 


creating a Zeeman effect in the absence of an anomalous magnetic moment term. The first two terms in the r.h.s. of 
(15) were found using the regularization procedure discussed in |3], and the vacuum term (14) was regularized with 
the help of Schwinger’s proper time scheme. 

The electron thermodynamic potential is 


n e = nr c + n™ ed = 


1 \eH\ 

4^/tt (2t r) 2 J_ c 


dp3 ^2 d(n) 


1/A2 S 


ds 
3/2 ' 


-,-sEi 


\eH\ 

47T 2 


/ OO 

dp:i {Pe -^e)C(Me ^e)|e=l? 

-oo 


(18) 

where the degeneracy factor d{n) = 2 — <5 n0 takes into account the lack of spin degeneracy of the LLL and the modes 
are given by the well-known spectrum of a free charged fermion in a magnetic field, 


E f . = e 


\Jrn 2 e +pl + 2\eH\n , 


e = ±. 


To find the expectation values of A f, bf, and (i/ we must solve the set of equations 




= 0, 




dAf ’ dbf 

together with the electrical neutrality condition 


= 0, 


d^ 2 E 


= 0, fe{u,d}, 


d^ 2 E 

d/J, e 


= 0. 


(19) 


( 20 ) 


( 21 ) 


The third equation in (20) is equivalent to the equation for the baryon density 3 ps = —dVl^E/dp. Its solution 
corresponds to a maximum of fl^ 2 E |23| . 

The contribution — Pfbf in (151 originates from the asymmetry of the spectrum at the LLL and is directly 
connected to the baryon charge anomaly [3]. Just as in the isospin-symmetric case the presence of this term 

favors nonzero values for bf for all /i > 0 within the region of validity of the model. Therefore, strictly speaking, 
for quark matter in the presence of a magnetic field, the formation of the CDW condensate is energetically favored 
from arbitrarily small to intermediate values of chemical potentials. This is clearly seen in the plots of Fig. [T] where 
the behavior of the magnitude and modulation of the condensates as functions of the baryon quark chemical are 
depicted for H = 2.5 x 10 18 G. The separation of the parameters bf,Af for each quark is a consequence of the 
neutrality condition that leads to different quark chemical potentials and hence different condensate solutions for 
different flavors. However, the behavior of the individual condensates with the chemical potential is similar. For small 
chemical potentials, the inhomogeneity is present, but the modulation is very small and the magnitudes of the up 
and down condensates coincide and are equal to the magnitude of the homogenous case. In contrast, in the region 
of interest for star applications, 330 MeV < p < 550 MeV, the two inhomogeneous condensates are quite distinct and 
robust. At very large chemical potential a competition may occur between the CDW solution and some form of color 
superconductivity, a topic worth being explored, but out of the scope of this paper. 
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B. Three-flavor model 

While the exact densities reached in the core of compact stars are unknown, it is likely that strange quarks might 
play a role in the thermodynamics of these systems. In particular, when building a hybrid EoS, if the transition 
from nuclear to quark matter occurs beyond the hyperon onset in the nuclear phase, a two-flavor description of quark 
matter will clearly lead to inconsistencies and the inclusion of a third flavor in our quark model might therefore be 
necessary. As in the two-flavor case, imposing the condition of /^-equilibrium leads to a split of the different flavor 
chemical potentials, which are given by 

2 1 1 

Mm — M — 2 Me ) Md — M T 2 Me > Ms — M A gMe ■ (22) 

We then consider a three-flavor version m of the neutral NJL model discussed above, with Lagrangian density 

A 3/) = Co + Ci + C v , (23) 


with 


A) = ^ (*7 + M7° + ip e (i^D^ + Me7° - ™e) ip e , 


(24) 


Ci = g s ^2 m a ^r + An 5 A a V0 2 


a =0 


(25) 


O 

C v = -G v Y [('</’7/A £ V) 2 + W’7 5 7A a V0 2 ] 


(26) 


a—0 


Here m = diag(m„, md , m s ) is the current mass matrix (for our calculations we again neglect the light quark masses, 
while choosing for the strange current mass a value of m s = 150 MeV), /t = diag(M miMcZjMs) is the flavor chemical 


potential matrix, A a are the Gell-Mann matrices in flavor space for a = 1, ... 8 , and Ao = 1. The covariant 

derivatives are defined as before, but with the replacement of the electric charge matrix by Q = diag{|e, — ^e, — ^e} 
for the quarks. When choosing our specific ansatz for the mean fields we follow [21] and allow only the light quark 
condensates to become spatially inhomogeneous with the same plane-wave form of Eq. ([8]), while implementing strange 
quarks as a homogeneous background of quasiparticles with constituent mass A s . Here again, the introduction of 
a nonzero Me breaks the isospin symmetry in the system, allowing for different values of the quark condensates of 
different flavors. As in the two-flavor case, we introduce quark number densities pj for each flavor, and replace each 
quark chemical potential by the effective one jlf = M/ — 4 GyPf- 

We have chosen not to include the six-fermion interaction generated by the instanton contribution in the three-flavor 
case [T§ because for isospin-asymmetric matter, with three condensates and two spatial modulations, it would give 
rise to coordinate-dependent terms in the mean-field Lagrangian due to the mixing of different condensates like for 
instance, ('4> u ' l Pu)(4 , d'4 , d)ips' t Ps, that would make our numerical calculations quite involved. Nevertheless, even though 
our simplified model is just a first attempt to tackle a very complicated problem, we do not expect that including 
the instanton term will change the physical picture significantly, as the mass-radius curves should be less sensitive to 
the degree of mixing than to the main features of the model, namely the asymmetry of the LLL Hamiltonian of the 
light quarks, and the existence of condensates with different magnitudes and modulations induced by the neutrality 
condition, all properties that are present whether there is mixing or not. 

Working in the mean-field approximation, we readily find that the spectrum of the u and d quarks is still given by 
Eq. (16 1, thus identical to the two-flavor case. The spectrum of the electrons is the same as before and the energies 


of the s quarks are given by 


E s = eyjA 2 s +pl + 2\e s H\n, e = ±. 
The thermodynamic potential for the three-flavor case is then given by 

n (3f) = q( 2/) + 7Vc ^ s + 


(A s -m s ) 2 (Ms — Ms) 2 


8G S 


8 G v 


(27) 


(28) 
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where fis given by Eq. (jT2J) and 

1 \e,H\ 


= 


4-A (2tt) 2 J_ c 


dp 3 ^2d{n) I 

n il/A 2 S 3/2 


-se: 


\e s H\ 
47r 2 


/ OO 

dp3(/U - E S )Q(J1 S - E s ) | e= i. 

-OO 


(29) 


The dynamical parameters of this model have to be determined from the equations 



dn( 3 E 



> 

II 

o 

dpf 

= 0, f = u,d,s, 

(30) 


dflW 




dbf 

-e 

s' 

II 

'•*-» 

o 

II 

(31) 






dp e 

= 0. 

(32) 


Since the s-quark condensate is homogeneous, its role is basically to decrease the number of electrons required for 
neutrality, but it does not produce any significant qualitative change in the solutions of the inhomogeneous condensates 
of the light quarks, which display the same behavior with the chemical potential as in the two-flavor case depicted in 
Fig. 0 


III. NUCLEAR MATTER IN A MAGNETIC FIELD 

For describing nuclear matter we employ the nonlinear Walecka model 12-51 . In the presence of an external magnetic 
field, it is characterized by the following Lagrangian [25] : 


E = ^2 Ei + ^2 + E m 

i B 


(33) 


where 


Ei = - mi)iln , (34) 

E b = ~ e B'iu - m B + 9<j B v ~ g^Bl^ ~ g pB r ■ pl^ipB , (35) 

E m = ^(d p ad^a - to 2 ct 2 ) -U(a) + ^m 2 ^^ - ^u liV w ,lv + ^m 2 p p p ■ p p - ^p P „p^ ■ (36) 


The sum is taken over baryons ( B ), considering the nuclear octect (protons, neutrons, and hyperons A, E _ , E + , E°, 
S - , and S°), and leptons ( l ), considering electrons and muons. The mesons (M) considered comprise the scalar er, 
isoscalar-vector uj p , and isovector-vector p p . They mediate the interactions between the baryon Dirac fields ipg. The 
lepton Dirac field is represented by ipi. Here, es is the electric charge of each baryon, t = (ti,t 2 ,T 3 ) denotes the 
isospin matrices, and mg is the baryon mass. The field tensors for the mesonic fields are given by u pl , = d p u) v — d v u} p 
and p^ = d p pl — d u pjj — g pB {pl l x pi), while U(a) = l/3bm n (g (7 NO') 3 + l/4c(g a Ncr) 4 is the scalar self-interactions, 
m n being the nucleon mass. 

In the mean-field approximation the mesonic fields <r, ujq , and pj] are assumed to acquire nonzero expectation values: 
(a) = a, (w 0 ) = w 0 , and (pj]) = pg. The mesonic masses are m a = 400 MeV, m u = 783 MeV, and m p = 770 MeV. 
The hyperon couplings to the mesonic fields are described as a fraction of that of the nucleons and are taken as 
x a H = 0.7 and x p h = x u h = 0.783 (xib = giB/gi , i = er,p,w). The parameters are chosen to reproduce a binding 
energy of —16.3 MeV and a symmetry energy coefficient of 32.5 MeV for saturated nuclear matter with compression 
modulus K = 300 MeV and effective baryon mass m* B = mg ~ gad = 0.7 ms- For this we adopt the following values: 
{ga/ma) 2 = 11-79 fm -2 , (p w /m w ) 2 = 7.149 fm -2 , (g P /m P ) 2 = 4.411 fm -2 , b = 0.002947, and c = —0.001070 (GM1 
parametrization). 

The thermodynamic potential at zero temperature is therefore: 


1 ( Pc 


n = — V fL — 

^ 2 V m u 

i=B,l 

Pi 3 y ( x P j i-M Pi i 

i=B 

pb — y ", x u ipi , 


Pb + l ( — 
2 V IXln 


1 1 

(gad) 2 + -bm n (g a a) 3 + -c(g a a) 4 


i= B 


1 

2 



(37) 

(38) 


(39) 








where the terms for charged particles (with dynamics modified by the filling of the Landau levels) and uncharged 
particles are given by 


a 


Pi 


neutral _ __ Tz 

3 (2?r ) 3 

^3 

neutral _ 


d A k 


_ k 2 ^ 

\J k 2 + to* 2 


37T 2 


n ?»r 9 e d = d{n) r dk 


n—0 

charged l e #1 vW . 

Pi = TTIT JL d ( n ) k fi 


Vk 2 


o*2 


27T 2 


n—0 


(40) 

(41) 

(42) 

(43) 


where 7 j is the degeneracy factor and the number density. The spin degeneracy of the Landau levels n is denoted 
as before by d{n) = 2 — S n Q. The sum is taken from the LLL to n maxi where n max is the nearest natural number 
equal to or less than [(ju 2 — m* 2 )/2\eH\] with pb = /iq — g u uJo — 9p T 3iPo denoting the effective chemical potential for the 
given fermion. If we write the effective mass of charged components as to * 2 = to * 2 + 2n|ei?|, the Fermi momentum 
of the particles becomes k 2 4 = p 2 — to* 2 . The baryon sum can be limited to protons and neutrons (GMln case) or 
include the full baryon octet (GMlnh case). The EoS for the two cases will start to differ at densities around 0.3-0.4 
fm ~ 3 or he ~ 1230 MeV, which mark the onset of hyperons |3 HI]. The thermodynamically favored values for the 
variational parameters in the model are obtained by solving the field equations describing the coupling of baryons to 
the mesons while considering chemical equilibrium and charge neutrality: 


( 5 ^) ^ 9<7<T ^ + 6m "^ CTCr ^ 2 + c ( fl, ' T<7 ) 3 


w 0 


Po 

i=B,l 

Pi 



E 9u 

TO 2 3tt 2 ’ 


9p k }B T 
^m 2 p 3n 2 3B ’ 

B P 


0, 


k^iPB T &iPe ) 


(44) 

(45) 

(46) 

(47) 

(48) 


where ej and Bi are the electric and baryonic charges of each component, associated with the corresponding electrical 
and baryonic chemical potentials, fi e and ps, respectively. 


IV. VARYING INNER MAGNETIC FIELD 

Most compact stellar objects are known to have very high values of surface magnetic fields, with white dwarfs in the 
range of 10 6 — 10 9 G, typical neutron stars 10 s — 10 12 G, and magnetars 10 14 — 10 15 G [2S]. While the magnetic field 
is definitely expected to rise in several orders of magnitude when going from the surface to the core, its actual inner 
profile is unknown. Estimates based on the virial theorem for stars made of quark matter give upper values of central 
fields of order 10 19 — 10 2O G 23- Other estimates found by solving the Einstein equations with axisymmetric and 
poloidal field configurations EH HO, or by applying the virial theorem to stars entirely composed of nuclear matter 
[32] . have led to the lower range 0.1 — 4.2 x 10 18 G. 

In all our derivations, we consider a static background magnetic field H pointing in the z-direction. This field 
influences the EoS both by altering the energy spectrum of the charged particles and by producing a splitting between 
the parallel and transverse pressures with respect to the field direction [29] ■ The pressure splitting is mostly due to 
the Maxwell term £em = — where F t , u = — d v A p , contribution to the Lagrangian. We shall work 

in a region of magnetic field strengths where the pressure splitting is ^ 10 %, so that the ambiguity associated with 
using the spherical TOV equations to obtain the mass-radius sequences is also small. Additionally, we neglect the 











9 



FIG. 2. Magnetic field profiles corresponding to Eq. (491 for different parametrizations (k, 7 ) in the chemical potential region 
under consideration. 


interaction of the magnetic field with the anomalous magnetic moment of the particles because its inclusion has been 
shown to have a negligible effect on the EoS of the dense medium [33] . 

In the stellar medium, the electric conductivity is very large; thus the magnetic flux is conserved. Hence, the 
magnetic field should be stronger as the density increases towards the core. Assuming a constant held magnitude 
throughout the star would then be a very crude approximation, which might introduce a significant bias in the 
resulting EoS. To avoid this, we consider a varying magnetic held in the star. A hrst implementation of a varying 
magnetic held inside the star was done in Ref. 14] and then used by several authors [34] . There, the neutron star 
was assumed to be composed entirely of nuclear matter and the magnetic held was expressed as a function of the 
baryon density, changing from a maximum central value to some lower surface value that was estimated from known 
magnetars. This ansatz, however, is not convenient to study hybrid stars on which the change from nuclear to quark 
matter occurs through a first-order transition. The jump that occurs in the baryon density due to the first-order 
transition would in turn produce an unphysical jump in the magnetic held. In such a case, a better way to mimic the 
inner varying held is to express it as a function of the baryon chemical potential, as proposed in m- Here we adopt 
the same approach and consider a medium-dependent value of H when calculating the EoSs of the nuclear and quark 
phases. With this aim, we employ the following ansatz: 


HM = H s + H c 


1 — e 


(49) 


where /iw = 938 MeV can be interpreted as the baryon chemical potential for nuclear matter at the crust of the star. 
The parameters k and 7 determine how fast the rise of H is with the baryon chemical potential from the surface to 
the core. Hg is the value for the surface magnetic field and He an estimate of the value at the core. The uncertainty 
in the knowledge of the actual inner profile of the magnetic field is expressed in (49) by the arbitrariness of the ( 7 , k) 
parametrization. While building hybrid stars, in Sec|~V)l we consider different values within acceptable ranges for 
these parameters, in order to determine the sensitivity of the EoS and the maximum stellar mass to their variation. 

In FigjSjwe show the magnetic field profiles for different sets of the (k, 7 ) parameters. Note that for the selected 
parametrization, the field decays quicker for densities corresponding to nuclear matter, while in the quark core it is 
almost constant. 


V. INHOMOGENEOUS QUARK MATTER IN THE CORE OF COMPACT STARS 

Since the first indications of the existence of 2 Mq stars [35] [36], there were claims [37, [38] that quark matter might 
have to be ruled out as a core phase. The problem was motivated by the fact that despite more than a decade of 
perturbative QCD predicting maximum stellar masses for quark stars larger than 2 Mq ESI SO], the expected densities 
of compact stars are not sufficiently high to validate a perturbative approach for the strong interaction. On the 
other hand, nonperturbative calculations based on simple QCD phenomenological models like NJL with four-fermion 
channels were found to produce too soft an EOS, incapable to stabilize a 2M 0 star against the gravitational collapse 
(see for example [41] and Fig. 8 in |15j ). However, when other interactions, also present in a dense medium of quarks 
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|42j , such as the diquark channel [53] and/or the vector channel, were taken into account, the EOS stiffened enough 
to meet the 2 Mq challenge. Because of these findings, quark matter was back in the competition as a possible core 
phase of massive stars (see for instance [551 - 155] '). One may wonder if this last conclusion could be challenged in turn 
by the presence of an inhomogeneous quark matter phase in the moderate density region. 

Exploring the suitability of the CDW phase for the core of neutron stars is particularly important, given that 
practically all neutron stars have nonzero magnetic fields, and once a magnetic field is present, the CDW phase 
is known to be energetically favored over the homogeneous chiral phases of quark matter, the chirally broken (at 
low p) and the chirally restored (at intermediate p). As discussed in Sec. II, the reason why the CDW solution is 
energetically favored is connected to the anomaly of the baryon charge produced by the spectrum asymmetry of the 
LLL. The charge anomaly leads to the term — in the thermodynamic potential that increases the pressure, 

rendering the EoS stiffer. Moreover, from the behavior of the parameters bf with the quark baryon chemical potential 
shown in Fig. [l] it is evident that this term will become more and more important with increasing p. What is unclear 
however is whether this term alone will be enough to compensate for the softening of the EoS typically associated 
to a transition to quark matter, or if other common stiffening factors such as vector interactions will also be needed. 
The only way to find out is through explicit calculations. Consequently, a main question we aim to explore in this 
section is the following: Can a magnetized hybrid star with a core of quark matter in the CDW phase sustain a star 
mass consistent with the 2 Mq observations for an acceptable range of the model parameters? 


A. Pressure splitting and Maxwell construction 

To find the transition point from nuclear to quark matter, we employ the Maxwell construction, prescribing that the 
transition between two phases 1 and 2 occurs at the same baryonic chemical potential, /i bi = /i B 2 , temperature, 
Ti = T 2 , and pressure, Pi = P 2 . The transition is then of first order and the density exhibits a discontinuity at the 
phase transition. It was shown in m that the Gibbs construction, for which the continuity of the electron chemical 
potential is also required, leads to very similar results for the macroscopic properties of a compact star, so that the 
choice of the Maxwell construction should be acceptable. 

The inclusion of a background magnetic field can in principle introduce a richer scenario in the construction of a 
hybrid EoS. Indeed, as mentioned in the previous section, when a system is subject to an uniform magnetic field along 
a specific direction, the pressure of the system develops a splitting in the directions parallel and perpendicular to the 
applied field. This splitting has to be taken into account in the EoS, which is then modified from the usual form into 

m 


P" = -n - ^, (50) 

tj2 

P ± = -n- HM + — , (51) 

i/2 

s = fl + fj,p + -y - , (52) 

where f 1 is the matter contribution to the thermodynamical potential evaluated at the physical minimum, p = —dfl/dp 
is the particle density, e the energy density and A4 = —dfl/dH the magnetization. 

In light of the pressure anisotropy, the Maxwell construction has to be generalized to require separately the equal¬ 
ity of the pressure components of the two phases. Nevertheless, taking into account that the leading term in the 
thermodynamic potential is ~ /i 4 , and that for nonferromagnetic media M < H, it follows that for the region under 
consideration, where H < p 2 , neglecting the magnetization energy (A iH) in the transverse pressure is a reasonable 
assumption. As a corroboration of these arguments, a direct calculation of the magnetization term for quark matter 
shows that for the range of p relevant for star applications 330 MeV< p < 500 MeV, HA4/Q never exceeds 4% (see 
Fig. §. while for nuclear matter HA4/H 2 does not exceed 3% [53] . Thus, the magnetization contribution can be 
neglected in the two phases. This, together with the fact that the magnetic field at the interphase is the same, reduces 
the Maxwell construction to 


^nuclear (/^tr) — f^quark(/^r) 5 


(53) 


with Unuciear given by Eq. (37) and f2 qua rk given either by Eq. (12) or Eq. (28), depending on whether the quark 
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p (MeV) 


FIG. 3. Magnetization contribution compared to the matter thermodynamic potential (Eq. as a function of the quark 

chemical potential for a two-flavor system in a fixed external magnetic field. The Haas-van Alphen oscillations are richer in the 
isospin-asymmetric matter here considered due to the difference in the maximum Landau levels reachable for the up and down 
quark. 


matter is composed of two or three flavors. Using this procedure, the transition chemical potential /i tr can be obtained 
from the solution of Eq. (53). 

Upon closer inspection, one can see that this criterion is however not entirely consistent: since the NJL model does 
not include gluonic degrees of freedom, its pressure will be completely blind to any effect related to confinement, 
while the nuclear model considered obviously deals exclusively with confined objects. In this sense, the transition 
occurring at /i tr should be interpreted as a “deconfinement” phase transition, which is not built from fundamental 
gluon dynamics, but only as an effective construction relying on the two phenomenological models involved. While 
a proper inclusion of confinement properties and the interplay between chiral and deconfinement phase transitions is 
clearly beyond the scope of this paper, in the spirit of previous works mm we attempt to incorporate these effects 
in a crude way through the introduction of a constant shift to the NJL vacuum pressure <5f2o, which will be treated 
as a free parameter. 

In the following, as done e.g. in |?7] for the case of homogeneous condensates, we consider two possible scenarios, 
one where strange quarks do not contribute to the thermodynamics of the star, and another in which they are included. 
For the first case (which we refer to as SU(2) case), we neglect hyperons in the nuclear EoS (GMln case) and consider 
only the two light flavors for the quark part. In order to have a consistent SU(2) description, the phase transition 
to quark matter must occur before the onset of hyperons. This limits the maximum value of the vector coupling in 
our calculations to Gy ~ 0.02Gs with = 0. For the SU(3) case, such a limitation is not present, although if 
Gy > 0.05Gs the transition chemical potential is so high that a quark matter core is never realized. The hybrid EoS 
for each case, using the parallel and perpendicular pressures, is shown in Fig. [4] We note that the phase transition 
happening at high chemical potentials would present a much more prominent density jump. It is also noticeable in 
the curve for GMlnh+SU(3) the softening of the EoS when the strange quark appears, as well as the stiffening of the 
quark EoS with increasing Gy. 


B. Masses and radii 

We now present numerical results for (M-R) sequences using the hybrid EoS obtained in the previous section, together 
with the Baym-Pethick-Sutherland [50) EoS for the star crust. As discussed above, the presence of strong magnetic 
fields breaks the star’s spherical symmetry and the usual TOV equations for obtaining the star’s structure are no 
longer valid m- However, one can consider a range of magnetic fields that is physically meaningful for compact stars 
and yet does not produce a sizable splitting in the pressure. An example of the pressure splitting profile inside the 
star is given in Fig. [5j As expected, such splitting of the pressures is more prominent with an increase in the field, 
which translates to higher densities given the density dependence of the magnetic held inside the star. One can see 
from Fig. [ 5 ] that if the held remains below 3 x 10 18 G the relative error associated with using one of the two pressures 
as representative of the whole interior of the star remains relatively small (< 10%). A similar scale has been found 
using the EoS of other models of dense quark matter m- Therefore, we work within the region of helds satisfying 
H < 3 x 10 18 G, so as not to invalidate the use of the spherical TOV equations. As a cross-check, we perform our 
calculations using both pressures, in order to make sure that the choice of one over the other does not lead to dramatic 
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FIG. 4. Equation of state (considering P± on the left panel and Py in the right one) for hybrid stars for different values 
of the vector repulsion as indicated. The nuclear EoS does not include hyperons when considering quark matter with only 
quarks up and down (GMln+SU(2) case) and includes hyperons when the strange quark is included in the quark phase as well 
(G Mlnh +SU(3) case). In the first case, the phase transition must happen at lower values of fiB- The magnetic field follows 
Eq. (491 and is taken with Hs = 1 x 10 15 G, He = 2.5 x 10 18 G, 7 = 2.5, and k = 12. The horizontal lines joining the jumps in 
the energy density produced by the first-order phase transitions are only graphical artifacts to connect the curves before and 
after the phase transitions. No shift in the vacuum value of the quark phase was introduced. 



FIG. 5. Parallel (lower curves) and perpendicular (upper curves) pressures for a two-flavor hybrid star as a function of the 
baryonic chemical potential. The surface magnetic field is taken as 1 x 10 15 G and the central field is 2.5 x 10 18 G for the full 
lines and 6.8 x 10 18 G for the dashed ones with 7 = 2.5 and k = 12. The curves end at the maximum value of the baryon 
chemical potential achieved in the interior of each configuration. 


changes in our results. 

Hybrid star M-R sequences with an inhomogeneous quark matter core are shown in Fig. [6] for values of surface 
magnetic field compatible with magnetars, Hs » lx 10 15 G and He ~ 2.5 x 10 18 G, using the perpendicular pressure. 
From the shape of these curves, the shrinking of the quark matter core with increasing Gy is clearly visible. The 
effective field in the center of the star is usually smaller than He as can be seen in Table [E] As expected, the inclusion 
of strangeness for both the nuclear and quark phases softens the EoS, so that the maximum mass is substantially 
reduced. In the case of GMlnh+SU(3) with Gy = 0.05, the phase transition happens for a very high value of the 
central density and the quark core turns out to be extremely small. 

As previously mentioned, we also investigated the effects of including a constant shift in the NJL model vacuum 
pressure, which we interpret as a contribution from confining effects (i.e. a bag constant). Following jlll 149j. this 
quantity is treated as a free parameter ranging between —17 < <5f2o < 0 MeV/fm 3 . It is clear from Fig. pUhat the 
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FIG. 6 . The mass-radius relations for a two-flavor (left) and three-flavor (right) inhomogeneous hybrid star considering Hs = 
1 x 10 15 G, He = 2.5 x 10 1 S G, 7 = 2.5, and k = 12. We employed the perpendicular pressure and the shifted vacuum pressure 
in units of MeV/fm 3 (see text for details) is indicated. The curves for stars composed entirely of nuclear matter in the GM1 
parametrization with no magnetic field are also shown for comparison. The mass constraints (±<r) of pulsars PSR J1614-2230 
and PSR J0348+0432 are shown as shaded regions. 



GMln (77=0) 

GMln+SU(2) 

GMlnh (77=0) 

GMlnh+SU(3) 

Gy/Gs 

- 

0 

0 

0 

0.02 

- 

0 

0.05 

7 

- 

2.00 

2.50 

3.00 

2.50 

- 

2.50 

2.50 

K 

- 

4.70 

11.70 

35.20 

11.70 

- 

11.70 

11.70 

H s (xlO 15 G) 

0 

1.00 

1-00 

1.00 

1.00 

0 

1.00 

1.00 

H c (xlO 18 G) 

0 

2.50 

2.50 

2.50 

2.50 

0 

2.50 

2.50 


2.39 

1.84 

1.85 

1.87 

1.96 

2.03 

1.70 

1.98 

Umax / He 

- 

0.78 

0.94 

1.00 

0.93 

- 

0.56 

0.92 


TABLE I. Maximum mass values in units of solar masses obtained using the perpendicular pressure for hybrid stars composed 
of nuclear matter (GM1) and two-flavor (GMln-\-SU{2)) and three-flavor (GMlnh + SU (3)) quark matter with inhomogeneous 
condensates for different values of the parameters determining the magnetic field profile inside the star and the quark repulsion 
strength ( Gy )• 


effect of the charge anomaly alone is not large enough to push the mass to the desired 2M 0 value. However, choosing 
a negative value of SQq results in a slightly stiffer EoS for quarks and pushes the phase transition from nuclear to 
quark matter to lower chemical potentials (albeit always larger than /r = 350 MeV). This allows the use of a larger 
value of the vector coupling, which will also help make the quark EoS stiffer. The left panel of Fig. [ 6 ]shows that with 
the combination of all these effects, for suitable values of Gy and <5Q 0 , it is possible to increase the maximum mass 
achieved in the two-flavor case. Although the mass increases only by a relatively small amount, it is enough to make 
it compatible with the observations of PSR J1614-2230 and PSR J0348+0432. 

In the GMlnh+SU(3) case, a lower value of the transition chemical potential and a stiffer quark phase could also 
be achieved by the use of a negative shift in the vacuum value (see the right panel of Fig. [ 6 ]). However, in this case, 
due to the small size of the quark core, the influence of these effects is greatly reduced compared to the two-flavor 
scenario. Of course, this is limited by the nuclear EoS used here: different models might allow for larger quark cores 
and make these effects more noticeable. 


Table [I] shows the maximum masses obtained using different parametrizations for the magnetic field in Eq. (491, 
as well as different values for the vector coupling. The results are obtained considering the perpendicular pressure. 
In this case, the maximum masses of the considered magnetized hybrid stars are smaller than those of pure nuclear 
matter at zero field (GMln (77=0)). 

In order to test the dependence of the maximum mass with the magnetic field profile inside the star, we have 
changed the values of 7 and n so as to have a faster (7 = 3.0 and k = 35.2) and a slower (7 = 2.0 and k = 4.7) 
increase of the field as one moves toward the center of the star (see Fig. [2|, while keeping the same values for He and 
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FIG. 7. Comparison of the mass-radius relations for a two-flavor hybrid star obtained when using the parallel and perpendicular 
pressures. The magnetic field is modeled with Hs = 1 x 10 1 S G, He = 2.5 x 10 18 G, 7 = 2.5, and k = 12. The values of Gy 
and shifted vacuum pressure in units of MeV/fm 3 are indicated. We note that the decrease in the maximum mass attainable 
is less than 5% when using the parallel pressure over the perpendicular one. 



GMln+SU(2) 

GMlnh+SU(3) 

Gy/Gs 

0 

0.02 

0.05 

0 

0.05 

SQ 0 (MeV/fm 3 ) 

0 

0 

-10 

0 

0 

Mmax ■ P .L 

1.85 

1.96 

2.05 

1.70 

1.98 

Mmax • -P|| 

1.77 

1.88 

1.95 

1.67 

1.89 


TABLE II. Maximum mass (in units of solar mass) allowed for hybrid stars composed of nuclear matter (GM1) and inhomo¬ 
geneous condensate obtained using the perpendicular and parallel pressures. The varying magnetic field profile is given by 
H s = 1 x 10 15 G, H c = 2.5 x 10 18 G, 7 = 2.5, and k= 12 . 


Hs for the case Gy — 0 and two-flavor matter. The difference in the corresponding maximum masses can be seen in 
Table [T] We note that increasing the rate of decay of the field toward the surface, the maximum mass increases very 
slightly. For the parametrization used, this only amounts to up to ~ 1% in comparison to the profile with 7 = 2.5 
and k = 11.7 used for all other calculations. Nevertheless, we point out that a more significant difference could be 
reached by allowing for a very steep increase/decrease of the magnetic field with /i. Whether such a profile would be 
compatible with the internal structure of compact stars is not known, so in the present work we restricted ourselves 
to the more conservative parametrizations shown in the table. We also point out that enforcing a larger value of 
He (although excessively strong central fields would lead to large star deformations and invalidate the use of the 
spherically symmetric TOV equations) is also expected to influence the maximum mass value in a more prominent 
way. 

The differences arising from the choice of the parallel or perpendicular pressures in the TOV equations are presented 
in Table [TT] and FigjTJ The difference in the maximum masses never exceeds 0.1 solar masses and the two curves differ 
from each other only at high densities, where the magnetic field strength is higher. Unfortunately, given the limitations 
of the method used here to obtain the star structure in light of the anisotropy introduced by the magnetic field, we 
cannot conclude which of these results could better represent the maximum mass for a highly magnetized compact 
object. We therefore simply treat the two results as upper and lower limits of an uncertainty band. In this context, 
it is interesting to point out that the inclusion of a magnetic field using an axisymmetrical geometry has been found 
to increase the star mass |31j . a result that might indicate that the actual physical situation in our study could be 
closer to the upper than the lower limit. Nevertheless, it should also be kept in mind that the conclusions of [3T] are 
not exempt of limitations either, as they were found considering a poloidal field configuration and disregarding the 
modification of the matter part of the EoS by the magnetic held. 

Taking all these elements into account, we find that with a relatively low value of Gy and a realistic value for He, 
well within the range where the pressure anisotropy is small enough to justify the application of the usual spherically 
symmetric TOV equations, it is possible to achieve maximum masses around 2 M 0 , a result compatible with the 
precise mass measurements of PSR J1614-2230 (M = 1.97 ± 0.04M o [5]) and PSR J0348+0432 (M = 2.01 ± 0.04M o 
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0 )- A quark matter core characterized by an inhomogeneous chiral condensate can thus be seen as a viable internal 
composition of these compact objects. 


VI. CONCLUDING REMARKS 

We studied the effects of the formation of inhomogeneous chiral symmetry breaking phases on the EoS of quark matter 
in a magnetic field, and the consequences for the masses of hybrid stars. To describe quark matter in the core, we 
considered a conventional NJL model with scalar and pseudoscalar four-fermion channels and included several elements 
relevant for astrophysical scenarios, such as electrical neutrality, /3 —equilibrium and vector repulsion. Considering 
a CDW ansatz for each light flavor, we saw that the spectral asymmetry of the LLL gives rise to a term in the 
thermodynamic potential that favors the formation of this inhomogeneous ground state and stiffens the pressure. 
Using this quark model together with the well-established nonlinear Walecka model, we built hybrid EoSs for stars 
with quark matter core in the CDW phase, and showed that such a configuration can support masses of ~ 2Mq. 
We investigated the sensitivity of these results on the parametrizations involved, and found that masses compatible 
with the recent PSR J1614-2230 and PSR J0348+0432 measurements can be obtained with reasonable values of the 
parameters chosen, although in the case where strangeness is allowed, the quark core would end up being significantly 
reduced in order to achieve higher masses. 

We find the fact that the realization of inhomogeneous quark matter in the core of compact stellar objects gives 
compatible results with the current mass observation a very encouraging outcome of our investigation, even when 
considering the still unsolved issue with the splitting of the pressures in the presence of a magnetic field. In light of 
this, it would definitely be interesting to investigate other effects of the formation of inhomogeneous condensates on 
the physics of compact stellar objects, such as transport and cooling properties, with the hope of finding stronger 
experimental signatures. 

One may wonder if the CDW condensate of the magnetized quark system could be washed out by fluctuations. Let 
us recall that in one-dimensional systems, the Mermin-Wagner theorem (51j forbids the existence of (1+1)-dimensional 
Nambu-Goldstone bosons, so that in such systems, condensate solutions that break chiral symmetry cannot be stable 
against quantum fluctuations about the condensate. In addition, single-modulated condensates in (3+l)-dimensional 
systems can be quite sensitive to temperature fluctuations because their pairing dynamics is often essentially one¬ 
dimensional. This possibility was recently investigated in [52], where it was found that the Nambu-Goldstone modes 
associated with the CDW condensate wash out the long-range order at finite temperatures, although they do support 
algebraically decaying long-range correlations, so the phase can still exhibit a quasi-one-dimensional order as in liquid 
crystals. 

While the analysis of our paper has not included the effects of fluctuations, we believe that there are good reasons 
to expect that in the case with magnetic field, the fluctuations will likely not be as effective in inducing the instability 
of the CDW condensate. The magnetic field is known to enhance chiral symmetry breaking through the dimensional 
reduction of the LLL fermions, the well-known mechanism of magnetic catalysis (see |53| for a recent review). The 
same dimensional reduction has proven to be essential to make the inhomogeneous CDW condensate energetically 
favored thanks to the asymmetry of the LLL spectrum in the inhomogeneous background. Now, the magnetic catalysis 
mechanism was recently questioned by claims that the dimensional reduction of the LLL fermions would translate into 
an effective dimensional reduction of the Nambu-Goldstone bosons, which in turn would hinder the stability of the 
chiral condensate, the so-called inverse magnetic catalysis [53]. But these claims were later challenged by the results 
of Ref. |55| that used a functional renormalization group approach to demonstrate that the constituent quark mass 
increases with the magnetic field at all temperatures and concluded that despite a strong anisotropy in the meson 
propagation, their fluctuations do not lead to the inverse magnetic catalysis claimed in [53]. In view of this, we expect 
that a similar behavior to the one found in [55J should occur in the case of the inhomogeneous chiral condensate. At 
this point however, we admit that ours are just hand-waving arguments that can be corroborated only by a thorough 
study of the fluctuations in the CDW phase in a magnetic field, an interesting task to be undertaken in the future. 

We remark that in order to obtain a first insight on the effect of inhomogeneous quark matter on the stellar 
EoS, several simplifying assumptions were made when building our model. One was to ignore the generation of the 
condensate associated with the anomalous magnetic moment of the chiral pairs. Such a condensate has been found in 
the presence of a magnetic field for the homogeneous background |56j and even in the case of color superconductivity 
in a magnetic field [57H59] , There is no reason to expect it is not present also in the inhomogeneous case. However, 
the magnitude of the magnetic moment condensate is typically small compared to the chiral condensate, except for 
extremely large magnetic fields, hence, as a first approximation it can be neglected. Another simplification was 
the omission of color-superconducting phases, which are expected to be the true ground state at asymptotically 
high densities. While we recall that the presence of a magnetic field has been shown to favor the formation of 
chiral crystalline phases, it is also known that a magnetic field leads to extra condensates in color superconductivity 
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[53 ESI, and makes the MCFL more stable than the regular CFL, so it is expected that a competition between color 
superconductivity and CDW may occur at intermediate densities in the presence of a magnetic field. Therefore, it 
is a pending and important task to explore whether the inhomogeneous phase can push or not the onset of color 
superconductivity in the presence of a magnetic held to densities higher than the ones considered in the present work. 
An explicit model calculation to address this question would be of course highly desirable. Additionally, an interesting 
question to consider is whether the incorporation of gluon effects through an extension of the NJL quark model to a 
gauged NJL quark model could lead to any sizable softening of the EoS of the system, as was recently found to occur 
in the case of CFL color superconductivity [60]. 

Finally, we recall that the inhomogeneous phase considered in this work is just one of the possible exotic phases that 
could be realized in dense matter. Another plausible candidate for the ground state of strongly interacting matter 
at intermediate densities, the so-called quarkyonic matter eh, is also characterized by a spatially varying chiral 
condensate |6‘Jl 1631 , although possibly with very different characteristics, particularly in a magnetic field background 
[64] , which could result in new unexpected effects on the EoS and on the transport properties of cold and dense quark 
matter. 
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